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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

© Verfahren zur Bestimmung der Dichte eines Gases unter Betriebsbedingungen 

(57) Die Dichte eines Gases unter Betriebsbedingungen 
kann mit den heut serienma&ig verfugbaren Betriebs- 
dichteaufnehmern lediglich mit einer Genauigkeit von 3% 
gemessen werden. Bei dem neuen Verfahren wird ein ei- 
nen Schwingkorper enthaltender Betriebsdichteaufneh- 
mer mit Hilfe eines Kalibriergases geeicht. Anschliefcend 
wird ein DichtemeSwert des Gases unter Betriebsbedi- 
nungen mit Hilfe des Betriebsdichteaufnehmers gemes- 
sen. Aufterdem wird die Schallgeschwindigkeit und die 
Viskositat sowohl des Gases als auch des Kalibriergases 
jeweils unter Betriebsbedingungen bestimmt. Die Dichte 
des Gases unter Betriebsbedingungen wird aus dem 
Dichte mefcwert mit Hilfe einer von der Schallgeschwin- 
digkeit und der Viskositat des Gases und des Kalibrierga- 
ses abhangigen Korrekturfunktion mit einem Fehler un- 
terhalb von 0,1% bestimmt. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren zur Bestimmung der Dichte eines Gases, insbesondere von Erdgas, unter Be- 
triebsbedingungen, wobei ein einen Schwingkorper enthaltenderBetriebsdichteaufnehmer mit Hilfe eines Kalibrieigases 
5 geeicht wird, und ein DichtmeBwert des Gases unter Betriebsbedingungen mit Hilfe des Betriebsdichteaufnehmers ge- 
messen wird. 

Als Schwingkorper konncn in den Betriebsdichteaufnchmcrn z. B. Schwinggabeln oder Schwingzy Under verwendet 
werden. Zur Messung eines DichtemeBwertes wird der Schwingkorper an seiner Oberflache dem zu untersuchenden Gas 
ausgesetzt. Bei Anregung des Schwingkorpers in dem BetriebsoUchteaufnehmer wird auch das umgebende Gas zu 
to Schwingungen angeregt. Es werden Druck- und Geschwindigkeitsschwankungen im Gas erzeugt, die sich mit Schallge- 
schwindigkeit fortpflanzen. Dadurch wird vom schwingenden System Energie an das umgebende Gas abgegeben. Die 
vom Schwingkorper abgestrahlte Energie wird im Gas in kinetische Energie umgewandelt. Aufgrund der von der Dichte 
des umgebenden Gases abhangigen Dampfung des Schwingkopers nimmt die Periodendauer der Schwingungen zu und 
deren Frequenz ab. 

15 Bei der Kalibrierung eines Betriebdichteaufnehmers mit einem Kalibriergas wird ein Zusammenhang zwischen der 
Dichte des Kalibriergases und der Frequenz des Schwingkorpers aufgestellt. 

Wird der Dichteauf nehmer anschlieBend zur Bestimmung der Dichte eines Gases bei Betriebsbedingungen verwendet, 
welches eine andere Zusammensetzung als das Kalibriergas aufweist, treten systematische MeBfehler bis zu ca. 3% auf- 
grund des Gasarteneinflusses auf. 

20 Aufgabe der vorliegenden Erfindung ist es, ein kostengiinstiges Verfahren zur Verfugung zu stellen, um die Dichte ei- 
nes Gases unter Betriebsbedingungen genauer bestimmen zu konnen. 

Diese Aufgabe wird bei einem gattungsgemaBen Verfahren erfindungsgemaB dadurch geldst, daB die Dichte oder 
Schallgeschwindigkeit sowohl des Gases als auch des Kalibriergases jeweils unter Normbedingungen erfaBt wird und die 
Dichte des Gases unter Betriebsbedingungen aus dem DichtemeBwert mit Hilfe einer von der Differenz der Dichten oder 

25 der Schallgeschwindigkeiten des Gases und des Kalibriergases unter Normbedingungen abhangigen Korrekturfunktion 
bestimmt wird. 

Eine auBerst einfache Korrekturfunktion lautet z. B. wie folgt: 

Pkorr = PBDA ' [1 + Kp n • Ap n ] (1), 

30 

wobei Kp^ eine Konstante oder eine Funktion von Pbda und/oder der Temperatur ist, und 

Ap u = Pn,mess — Postal- 

35 Genauso einfach kann die Dichte des Gases mit Hilfe der folgenden Korrekturfunktion bestimmt werden: 
Pkorr = Pbda • [1 + K^ • Ac n ] (2), 

wobei Kc n eine Konstante oder eine Funktion von Pbda und/oder der Temperatur ist, und 

40 

Ac n = C n mess — C n jkaj. 

Eine Weiterbildung der Erfindung, die eine noch etwas genauere Bestimmung der Dichte des Gases ermoglicht, ist da- 
durch gekennzeichnet, daB auBerdem der Brennwert des Gases und des Kalibriergases unter Normbedingungen erfaBt 
45 wird und daB die Korrekturfunktion, die zur Bestimmung der unter Betriebsbedingungen gemessenen Dichte verwendet 
wird, vom Brennwert des Gases und des Kalibriergases unter Normbedingungen abhangig ist. 
Beispielsweise kann die Korrekturfunktion in diesem Fall die folgende Form haben: 

Pkon- = Pbda • [1 + Kp n • Ap n + K Hn • AH n ] (3), 

50 

wobei Kp n und Kh„ Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, und 

Apn = pn,raess ~ Pnjsal Und 
AH n = H n mess — Hnjcaj. 

55 

Stehen bei einer speziellen Anwendung bereits Schallgeschwindigkeitswerte fur das Gas und das Kalibriergas unter 
Normbedingungen zur Verfugung, so kann die Dichte des Gases alternativ mit Hilfe der folgenden Korrekturfunktion be- 
stimmt werden: 

60 p k orr = PfiDA * [ 1 + * Ac n + K Hn ■ AH n ] (4), 

wobei Kc n und K Hn Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, und 

Ac n = c n>raess - c nJcal und 
65 AH n = H nracs5 - H^kaj. 

Die MeBgenauigkeit zur Bestimmung der Dichte des Gases kann dadurch weiter erhoht werden, daB auBerdem die 
Temperatur des Gases unter Betriebsbedingungen (T raess ) und die Temperatur des Kalibriergases beim Kalibrieren (Tfcai) 
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erfaBt werden und daB die Korrekturfunktion, die zur Bestimmung der unter Betriebsbedingungen gemessenen Dichte 
verwendet wird, von der Betriebstemperatur des Gases und der Kalibriertemperatur des Kalibriergases abhangig ist. 

NaturgemaB laBt sich dabei die hochste Genauigkeit dadurch erzielen, daB sowohl der Brennwert des Gases und des 
Kalibriergases unter Normbedingungen als auch die Betriebstemperatur und die Kalibriertemperatur erfaBt werden. 

Es hat sich gezeigt, daB die Dichte des Gases besonders zuverlassig mit Hilfe der folgenden Korrekturfunktion be- 5 
stimmt werden kann: 

Pkorr = Pbda • [1 + Kp n • Ap n + K„ n • AH n ] + Kt • AT (5), 

wobei Kp n und Kh q Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, Kt eine Konstante oder eine 10 
Funktion von Pbda isU und 

Ap n = Poorness ~ Pnjcal, 

AH n = H n mess - H Il kil ] und 

AT = T raess - Tfcai. 15 

Genauso zuverlassige Ergebnisse fur die Dichte konnen mit Hilfe der folgenden Korrekturfunktion bestimmt werden: 
Pkorr = Pbda ■ [1 + • Ac n + K Hn • AH n ] + K T • AT (6), 



wobei K^ und K Hn Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, K T eine Konstante oder eine 
Funktion von Pbda ist, und 
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Ac n — Cn mess ~ c njcal, 

AH n = H n>mcss - H nJcai und 25 
AT = T raess - Tkai. 

Eine Weiterbildung der Erfindung, mit der die Dichte des Gases unter Betriebsbedingungen noch genauer bestimmt 
werden kann ist dadurch gekennzeichnet, daB auBerdem der Kohlenstoffdioxidanteil des Gases und des Kalibriergases 
( x co2,mess» *C02,kai) erfaBt wird und daB die Korrekturfunktion, die zur Bestimmung der unter Betriebsbedingungen ge- 30 
messenen Dichte verwendet wird, von dem Kohlenstoffdioxidanteil des Gases und des Kalibriergases abhangig ist. 

Bei einem bevorzugten Ausfuhrungsbeispiel wird die Dichte des Gases dabei mit Hilfe der folgenden Korrekturfunk- 
tion bestimmt: 

Pkorr = Pbda • [1 + Kp n • Ap n + K Hn • AH n + Kco2 • Ax C0 2] + K T • AT (7), 35 

wobei Kp n , Kh„ und Kco2 Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, K T eine Konstante oder 
eine Funktion von Pbda un d 

Apn = p n ,racss ~ Pnjcal, 40 
ss ~ H^jjaj, 

- Tka, und 

AXC02 = Xc02,mess ~ x C02 T kal. 

Alternativ kann die Dichte des Gases mit der folgenden Korrekturfunktion bestimmt werden: 45 
Pkorr = Pbda • [1 + Kc a ■ Ac n + K Hn • AH„ + Kco2 ■ Axco 2 ] + K T ■ AT (8), 

wobei Kc n , Kh d und Kco2 Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, Kt eine Konstante oder 
eine Funktion von Pbda i sl > un d 50 

Ac n = c niTiess — Cn^ai, 
AH n — H n raess H n> j^ai, 
AT = T raess - T kal und 

AXC02 = Xc02,raess ~ XC02Jcai- 55 

Vorzugsweise wird die Variable K T mit Hilfe der folgenden einfachen Funktion bestimmt: Kr = Kt 0 + Pbda ■ Kt 17 wo- 
bei K To und K Tl Konstanten sind. 

Die oben genannte Aufgabe der vorliegenden Erfindung wird ferner von einem gattungsgemaBen Verfahren gelost, das 
dadurch gekennzeichnet ist, daB die Schallgeschwindigkeit und die Viskositat sowohl des Gases als auch des Kalibrier- 60 
gases jeweils unter Betriebsbedingungen (c mcss , c^, T| mcss , T]kai) bestimmt wird und die Dichte des Gases unter Betriebs- 
bedingungen aus dem DichtemeBwert mit Hilfe einer von der Schallgeschwindigkeit und der Viskositat des Gases und 
des Kalibriergases abhangigen Korrekturfunktion bestimmt wird. 

Die Schallgeschwindigkeit und die Viskositat des Gases und des Kalibriergases konnen entweder unrnittelbar gemes- 
sen werden oder aus anderen bekannten MeBgroBen abgeleitet werden. Die letztgenannte Moglichkeit wird in der Praxis 65 
aus Kostengriinden bevorzugt. 

Ein bevorzugtes Ausfuhrungsbeispiel der Erfindung ist dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe 
der folgenden Korrekturfunktion bestimmt wird: 
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20 



25 



1 + 



Kf 



+ K„ 



Pkorr f^3DA 



JfCal. 



1 + 



Kt 



i p.. 



ftcorr 



[(Kto + Pbda • K n ) • ( Tmes8 • T^J] , 



(9) 



10 wobei x: Periodendauer des Schwingkorpers, 
T: Temperatur und 
Kc, Kq, Kto, Kti Konstanten sind. 

Die mit Hilfe dieser Korrekturfunktion erreichbare Zuverlassigkeit bei der Bestimmung der Dichte des Gases unter 
Betriebsbedingungen ist betrachtlich. Dabei hat es sich bewahrt, zur Bestimmung der schallgeschwindigkeitsabhangigen 
15 Konstante Kc 2 Ethan als Referenzgas zu verwenden. 

Fur manche Anwendungsfalle, z. B. bei der Dichtebestimmung von Erdgas, kann eine entsprechend hohe Genauigkeit 
mit Hilfe der folgenden vereinfachten Korrekturfunktion erreicht werden: 



Pkorr = Pbi 



DA 



1 + 



K 



"Kal 



(■ 



+ K„ 



V Pbda 




Kal J J 



+ [fa™ + Pbda • k t1 ) • (i; e33 - TfCal )] 



(10) 



wobei X: Periodendauer des Schwingkorpers, 

T: Temperatur und 
30 Kc, Kq, Kjo, K T1 Konstanten sind. 

Grundsatzlich ist es vorteilhaft, den Betriebsdichteaufnehmer mit einem Kalibriergas zu kalibrieren, welches ungefahr 

die gleiche Zusammensetzung aufweist, wie das MeBgas, dessen Dichte anschlieBend gemessen werden soil. Soli die 

Dichte von Erdgas bestimmt werden, ist es vorteilhaft, daB Methan als Kalibriergas verwendet wird. Die Dichte von Me- 

than ist mit sehr hoher Genauigkeit im Stand der Technik bekannt. AuBerdem besteht Erdgas in der Regel zu wenigstens 
35 80% aus Methan. Das Erdgas aus den GUS-Staaten enthalt sogar mehr als 98% Methan. Somit bildet eine Kalibrierung 

mit Methan eine gute Basis fur eine genaue Dichtebestimmung von Erdgas. 

Vorzugsweise werden Kc 2 und K n durch Kalibrierung mit einem zweiten Kalibriergas neben Methan, z. B. Argon, 

Ethen oder Stickstoff, bestimmt und werden als typenspezifische Konstanten behandelt, welche lediglich vom Betriebs- 

dichteaufnehmertyp abhangen. 

40 Es hat sich bewahrt, daB die Schallgeschwindigkeit des Gases und/oder des Kalibriergases als Funktion der Dichte un- 
ter Betriebsbedingungen und der Dichte unter Normbedingungen und/oder der Temperatur und/oder des Brennwertes 
und/oder des Kohlenstoffdioxidanteils bestimmt wird. Auch hier gilt das Prinzip, daB die Schallgeschwindigkeit um so 
genauer bestimmt werden kann, je mehr MeBgroBen in die Funktion zur Bestimmung der Schallgeschwindigkeit einge- 
hen. Die Genauigkeit bei der nachfolgenden Bestimmung der Dichte des Gases unter Betriebsbedingungen steigt dem- 

45 entsprechend. 

Eine Weiterbildung der Erfindung ist dadurch gekennzeichnet, daB bei Bestimmung der Dichte von Erdgas als Kali- 
briergas Methan verwendet wird, daB die Schallgeschwindigkeit von Methan (CCH4) nach folgender Gleichung berechnet 
wird: 

50 c CH 4 = A 0 + (A! + A 2 . T raess ) • T raess + (A 3 + A4 • Pbda + A 5 • p 2 BDA ) • Pb da' + (A* + A 7 • T raess ) • T raess • p BDA (11) 

und die Schallgeschwindigkeit des Erdgases (c racss ) nach folgender Gleichung bestimmt wird: 

C raess = Cch4 + [(Bl + B 2 • Ap n + B 3 • T raess ) + (B 4 + B 5 • T mess + B 6 • Pbda + B 7 • Pbda ■ ApJ • Pbda] * A n 
55 + [(Ci + C 2 • T raess ) + (C 3 + C4 • Pbda + C 5 ■ T racss ) • Pbda] • AH D + [D t • Pbda] • X C o 2 (12), 

wobei A0-A7, B1-B7, C1-C5 und Di Konstanten sind, Xco2 den Kohlenstoffdioxidanteil im Erdgas angibt, 

APn = Pn,mcs5 " Pn.OL, Und 

60 AH n = Harness - H n ,CH4. 

Die Gleichungen (11) und (12) sind sehr einfach aufgebaut und hinreichend genau. Mit ihnen kann die Schallge- 
schwindigkeit mit einem Fehler unterhalb von 0,1% bestimmt werden. 

Vorteilhafterweise wird die Viskositat des Gases und/oder des Kalibriergases als Funktion der Dichte unter Betriebs- 
65 bedingungen und der Dichte unter Normbedingungen und/oder der Temperatur und/oder des Brennwertes bestimmt. Ge- 
hen alle MeBgroBen in die Funktion zur Bestimmung der Viskositat des Gases bzw. des Kalibriergases ein, laBt sich die 
Viskositat des Gases mit der hochsten Genauigkeit bestimmen. Die Genauigkeit bei der Bestimmung der Dichte des Ga- 
ses ist entsprechend hoch. 



4 



DE 198 23 193 A 1 

Beispielsweise kann, wenn bei Bestimmung der Dichte von Erdgas als Kalibriergas Methan verwendet wird, die Vis- 
kositat von Methan verwendet wird, die Viskositat von Methan OICH4) nach folgender Gleichung berechnet werden: 

TICH4 = Ao + A t • T + (A 2 + A 3 • Pbda) • Pbda (13) 

5 

und die Viskositat des Erdgases (T| mC ss) nach folgender Gleichung besdmmt werden: 

Tlmess = T\CH A + Bi (1 — Pbda) ' Ap n + B 2 ■ AH n (14) 

wobei A0-A3, Bi und B 2 Konstanten sind, 10 

Ap n = p n ,racss ~ Pn,CH4 und 

AH n = H n>mess - H n ,CH 4 - 

Auch die Gleichungen (13) und (14) sind sehr einfach aufgebaut. Die mit ihnen bestimmten Viskositaten sind nur mit 15 
einem Fehler deutlich unterhalb von 5% behaftet und somit hinreichend genau fur die Korrektur des DichtemeBwertes. 

Es hat sich bewahrt, daB zurn Auffinden der Parameter und Konstanten geeigneter Funktionen, insbesondere der Kor- 
rekturfunktionen, mehrere MeBzyklen mit Referenzgasen durchgefuhrt werden, deren Zusammensetzungen bekannt 
sind. Als Referenzgleichung zur Berechnung der Schallgeschwindigkeit kann z. B. die AGA-8-DC92-Zustandsglei- 
chung verwendet werden. Voraussetzung fur deren Anwendung ist jedoch die exakte Kenntnis der Zusammensetzung der 20 
Referenzgase. Als Referenzgleichung zur Berechnung der Viskositat von Erdgasen kann z. B. der Ansatz von Herning 
und Zipperer gewahlt werden. Vorteilhafterweise sollten die Referenzgase den Bereich der Zusammensetzungen abdek- 
ken, in denen die Zusammensetzungen der spateren MeBgase voraussichtlich liegen. 

Im folgenden wird die Erfindung anhand mehrerer in den Zeichnungen dargestellten Ausfuhrungsbeispielen naher er- 
lautert. In der Zeichnung zeigen: 25 

Fig, 1 acht Diagramme, in denen der DichtemeBfehler gegen die Dichte ftir sechs reine Gase und vier verschiedene 
Betriebsdichteaufnehmer bei T = 10°C aufgetragen ist, und zwar in der linken Spalte bei Anwendung des bekannten Ver- 
fahrens und in der rechten Spalte bei Anwendung des erfindungsgernaBen Verfahrens; 

Fig. 2 acht Diagramme, in denen der DichtemeBfehler gegen die Dichte fur sechs Erdgase und die vier verschiedenen 
Betriebsdichteaufnehmer bei T = 10°C aufgetragen ist, und zwar in der linken Spalte bei Anwendung des bekannten Ver- 30 
fahrens und in der rechten Spalte bei Anwendung des erfindungsgernaBen Verfahrens; 

Fig. 3 acht Diagramme, in denen der DichtemeBfehler fur Ekofisk-Erdgas und die vier verschiedenen Betriebsdichte- 
aufnehmer bei verschiedenen Temperaturen gegen die Dichte aufgetragen ist, und zwar in der linken Spalte bei Anwen- 
dung des bekannten Verfahrens und in der rechten Spalte bei Anwendung eines zweiten Ausfuhrungsbeispiels des erfin- 
dungsgernaBen Verfahrens; 35 

Fig. 4 sechs Diagramme, die die Wirkungsweise des erfindungsgernaBen Verfahrens gemaB Fig. 1 und 2 veranschau- 
lichen; 

Fig. 5 sechs Diagramme, die die Wirkungsweise eines dritten Ausfuhrungsbeispiels des erfindungsgernaBen Verfah- ■■ 
rens veranschaulichen. 

In Fig. 1 ist der DichtemeBfehler gegen die Dichte fur sechs reine Gase bei T = 10°C aufgetragen, und zwar unterein- 40 
ander fur vier verschiedene Typen von Betriebsdichteaufnehmern. In der linken Spalte sind die DichtemeBfehler bei der 
Dichtemessung gemaB dem Stand der Technik und in der rechten Spalte die verbleibenden Restfehler bei Einsatz des er- 
findungsgernaBen Verfahrens dargestellt. Bei dem dargestellten Ausfuhrungsbei spiel wurde die Dichte unter Verwen- 
dung der folgenden Korrekturfunkdon besdmmt: 



Pkorr = Pi 



BDA 



1 + 



^Kal ) 

(15) 



45 



50 



Die Betriebsdichteaufnehmer wurden vor der Messung bei 10°C mit Methan kalibriert. Es wurden anschlieBend die 
reinen Gase (N2, CO2, C2H6, C2H4, Ar, Ne) untersucht, da sie jeweils charakteristische Eigenschaften aufweisen und ei- 
nen groBen Schallgeschwindigkeits- und Viskositatsbereich abdecken. Wie deudich sichtbar ist, gelingt es mit dem er- 
findungsgernaBen Verfahren den systematischen DichtemeBfehler von bis zu 3% auf weniger als 0,1% zu reduzieren. Mit 55 
dem neuen Verfahren kann der DichtemeBfehler somit wirksam auf ungefahr 1/30 des urspriinglichen MeBfehlers ge- 
senkt werden. 

Fig. 2 unterscheidet sich lediglich dadurch von der Fig. 1, daB in den Diagrammen die MeBwerte fur sechs Erdgase bei 
T = 10°C aufgetragen sind. Die Betriebsdichteaufnehmer wurden wie in Fig. 1 mit Methan bei 10°C kalibriert. Wie den 
Diagrammen klar zu entnehmcn ist, kann der systemadsche MeBfehler bei Erdgasen von bis zu ca. +0,5% aus weniger 60 
als 0,05% reduziert werden, also auf weniger als 1/10 des urspriinglichen MeBfehlers. Das erfindung sgemaBe Verfahren 
ist somit auch fur Erdgase sehr wirksam. 

Die Betriebsdichteaufnehmer werden in der Regel nur bei einer konstanten Temperatur, z. B. 10°G, kalibriert. Weicht 
die Betriebstemperatur von der Kalibriertemperatur ab, treffen folglich durch den Temperaturunterschied bedingte zu- 
satzliche DichtemeBfehler auf. Mit dem in Fig. 3 veranschaulichten Ausfuhrungsbei spiel der Erfindung konnen auch 65 
diese zusatzlichen DichtemeBfehler behoben werden. In Fig. 3 sind auf der linken Seite die systemadschen DichtemeB- 
fehler fur Ekofisk-Erdgas bei den drei Temperaturen, T = 0°C, 10°C und 20°C fur vier verschiedene TVpe n von Betriebs- 
dichteaufnehmer aufgetragen. In den Diagrammen auf der rechten Seite sind die Restfehler bei Anwendung des erfin- 
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dungsgemaBen Verfahrens eingezeichnet. Bei diesem Ausfuhrungsbeispiel wurde die Korrekturfunktion (10) verwendet: 



Pkorr — Pbda 



1 + 



K 




Kal /J 



10 



L5 



20 



25 



30 



35 



40 



45 



50 



55 



60 



65 



+ [(K TO + Pbda ' K n ) • (^ ess - T Kal )] , 



wobei t: Periodendauer des Schwingkorpers, 

T: Temperatur und 

Kc, K^, Kto, K T1 Konstanten sind. 

Wie zu sehen ist, konnen auch in diesem Fall die DichtemeBfehler von rund +0,5% auf weniger als 0,05% reduziert 
werden. 

In Fig. 4 ist die Wirkungsweise der Ausfuhrungsbeispiele gemaB den Fig. 1 und 2 veranschaulicht. 
Die Korrekturfunktion (9) kann einfach in die folgende Form gebracht werden: 



BDA 



- Pi 



korr 



Pkor 



= Kt 



+ K. 



(16) 



wobei c* ein Schallgeschwindigkeitsparameter ist, namlich 

X 2 lc 2 D C 2 J ' 

^mess Kkorr Heal' 



(17) 



und TJ* ein Viskositatsparameter ist, namlich 



f airness 
I V Pkorr 



BDA 



Pko 



V Pbda> 



(18) 



Hierbei sind K c 2 und K n typenspezirische Aufnehmerkonstanten. 

Im Diagramm 1 von Fig, 4 ist der DichtemeBfehler gegen die Dichte fur sechs reine Gase aufgetragen. Diagramm 2 
unterscheidet sich dadurch von Diagramm 1, daB der DichtemeBfehler statt gegen die Dichte gegen den Schallgeschwin- 
digkeitsparameter c* aufgetragen ist. Wie dieses Diagramm bestatigt hangt der DichtemeBfehler linear von dem Schall- 
geschwindigkeitsparameter c* ab. Die eingezeichnete Gerade wurde mit Hilfe der verschiedenen MeBwerte angepaBt 
und liefert als Steigung die typenspezirische Aufnehmerkonstante Kc 2 . 

Im Diagramm 3 der Fig. 4 ist der Restfehler der Dichte aufgetragen, der sich ergibt, wenn die Korrektur beziiglich des 
Schallgeschwindigkeitsparameters bereits durchgefuhrt wurde. Dieser Restfehler ist gegen den Schallgeschwindigkeits- 
parameter c* aufgetragen. Wahrend die Werte fur Ethylen, Argon und Ethan bereits nur noch geringe Restfehler aufwei- 
sen, sind die Dichtewerte von StickstofF, Kohlenstoffdioxid und Neon mit noch relativ hohen Fehlern behaftet. Wie den 
Diagrammen 4 bis 6 zu entnehmen ist, konnen diese Restfehler mit Hilfe des Viskositatsparameters T|* deutlich reduziert 
werden. 

Im Diagramm 4 der Fig. 4 ist der Restfehler gegen den Viskositatsparameter T|* aufgetragen. Wie aufgrund der Funk- 
tionsgleichung zu erwarten ist, lassen sich die Restfehler in diesem Diagramm ebenfalls gut durch eine Gerade approxi- 
mieren. Die Steigung der Geraden liefert die typenspezirische Aufnehmerkonstante Kq. 

Im Diagramm 5 ist der Restfehler der Dichte nach Korrektur des Viskositatseffekts mit Hilfe des 2. Terms von Glei- 
chung (16) gegen den Viskositatsparameter T|* aufgetragen. Nachdem sowohl der Schallgeschwindigkeitseffekt als auch 
der ViskositatsefTekt korrigiert wurden, sind die Dichlen der sechs reinen Gase lediglich noch mit Restfehlem unterhalb 
von 0,1% behaftet. Dieser auBerst geringe Restfehler der Dichte ist im Diagramm 6 der Fig. 4 schlieBlich erneut gegen 
die Dichte aufgetragen, um einen besseren Vergleich mit dem DichtemeBfehler gemaB Diagramm 1 zu ermoglichen. 

Fur Erdgase kann die Dichte gemaB dem folgenden vereinfachten Ausfuhrungsbeispiel der Erfindung bestimmt wer- 
den. Es wurde festgestellt, daB die Schallgeschwindigkeit von Erdgas und die Viskositat von Erdgas uberraschend gut 
mit den bciden folgenden Gleichungen bestimmt werden kann: 

c E = cch4 (T,p) + Ac E (T,p,p n ) + Ac E (T,p,H n ) + Ac E (p,x C02 ) (19), 

TIE = T|CH4 (Xp) + ATtefoPn) + Al\ (Hn) (20), 

wobei die Ac E und At| e lediglich von den in den Klammern angegebenen MeBgroBen abhangige Korrekturterme bezeich- 



Werden diese Nahrungsterme fur die Schallgeschwindigkeit und die Viskositat des Erdgases in die Gleichung (10) ein- 
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gesetzt, ergibt dies naherungsweise die Funktion gemaB Gleichung (8): 
Pkorr = Pbda • [1 + Kp n • Ap n + K Hn - AH n + Kcoz • Axcozl + K T • AT, 

wobei K^, Kh d und Kco2 Konstanten, Funktionen von Pbda und/oder der Temperatur sind, 5 
Kt eine Konstante oder eine Funktion von Pbda ist, 

Apn = pn,racss — pnjcal' 
AH n = H ntraess - Hn^jcai, 

AT = T mcss - Tfcai und Ax C 02 = XC02,mess - x C02,kal- 10 

Bei dem in Fig. 5 dargestellten besonders einfachen Ausfiihrungsbeispiel der Erfindung ist Kco2 = 0 und Kp n sowie 
K Hn sind Konstanten und fiir K T gilt: K T = K To + Pbda ' K Tl , wobei Kt 0 und K Tl Konstanten sind. 

In Fig. 5 ist die Wirkungsweise dieses vereinfachten Verfahrens zur Bestimmung der Dichte von Erdgas fur sechs Erd- 
gase dargestellt. Im Diagramm 1 ist der DichtemeBfehler gegen die Dichte ohne Anwendung des erfindungsgemaBen 15 
Verfahrens aufgetragen. Wie dem Diagramm zu entnehmen ist, ist die Dichte mit einem Fehler von bis zu 0,5% behaftet. 

Im Diagramm 2 der Fig. 5 ist der DichtemeBfehler gegen die Dichte des Erdgases bei Normbedingungen aufgetragen. 
Wie aufgrund der Korrekturfunktion zu erwarten ist, lassen sich die MeBwerte in dieser Darstellung gut durch eine Ge- 
rade annahem. Die Steigung der Geraden liefert die Konstante Kp n . 

Im Diagramm 3 ist der Restfehler der Dichte nach Korrektur der Dichtewerte mit Hilfe des zweiten Terms der Glei- 20 
chung (8) dargestellt. Der Restfehler betragt fur alle sechs untersuchten Erdgase nur noch weniger als 0,1%. Im Dia- 
gramm 4 ist dieser Restfehler der Dichte gegen den Brennwert des jeweiligen Erdgases unter Normbedingungen darge- 
stellt. Auch hierbei lassen sich die bezuglich der Normdichte bereits korrigierten MeBwerte gut durch eine Gerade ap- 
proximieren. An der Steigung der Geraden laBt sich die Konstante K Hn ablesen. 

Im Diagramm 5 ist der Restfehler der Dichte nach Korrektur mit Hilfe des dritten Terms gemaB Gleichung (8) darge- 25 
stellt. Der verbleibende Restfehler betragt nur noch weniger als 0,05%. 

Im Diagramm 6 ist dieser Restfehler des erfindungsgemaBen Verfahrens gegen die Dichte dargestellt, um einen besse- 
ren Vergleich mit dem bisherigen MeBfehler gemaB Diagramm 1 zu errndglichen. 

Eine Korrektur bezuglich der Temperatur (4. Term der Gleichung 8) war in diesem Fall nicht erforderlich, da die MeB- 
temperatur mit der Kalibriertemperatur ubereinstimmte. 30 

Mit Hilfe des erfindungsgemaBen Verfahrens kann die Dichte je nach Anwendungsfall mit der gewiinschten Genauig- 
keit bestimmt werden. Die zur Korrektur des DichtemeBwertes benotigten Konstanten und Parameter der Korrekturfunk- 
tionen miissen fiir jeden Betriebsdichteaufnehmertyp lediglich einmal mit Hilfe von Referenzgasen bekannter Zusam- 
mensetzung bestimmt werden. 

Im Rahmen der Erfindung sind zahlreiche Weiterbildungen moglich. Beispielsweise kann auch Sticks toff als Kali- 35 
briergas verwendet werden. Die Kalibrierung sollte moglichst in der Mitte des iib lichen Betriebsternperaturbereiches lie- 
gen, jedoch konnen auch andere Kalibriertemperaturen gewahlt werden. Die zur Bestimmung des DichtemeBfehlers be- 
notigten exakten MeBwerte fiir die Dichte konnen bei bekannter Zusammensetzung entweder berechnet werden oder mit 
einem entsprechend kostenaufwendigen hochprazisen Betriebsdichteaufnehmer gemessen werden. 



Patentanspriiche 



3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Korrek- 
turfunktion bestimmt wird: 

Pkon- = Pbda • U + Kc n • AcJ, 

wobei Kc a eine Konstante oder eine Funktion von Pbda und/oder der Temperatur ist, und 

Ac n = C n raess — Cnjcaj. 



40 



1. Verfahren zur Bestimmung der Dichte eines Gases, insbesondere von Erdgas, unter Betriebsbedingungen, wobei 
ein einen Schwingkorper enthaltender Betriebsdichteaufnehmer mit Hilfe eines Kalibriergases geeicht wird, und ein 
DichtemeBwert des Gases unter Betriebsbedingungen (Pbda) nut Hilfe des Betriebsdichteaufnehmers gemessen 45 
wird, dadurch gekennzeichnet, daB 

die Dichte oder Schallgeschwindigkeit sowohi des Gases als auch des Kalibriergases jeweils unter Normbedingun- 
gen (p n ,raess> Pnjoi. c n>ra ess, c nJtal ) erfaBt wird und 

die Dichte des Gases unter Betriebsbedingungen aus dem DichtemeBwert mit Hilfe einer von der Differenz der 
Dichten oder der Schallgeschwindigkeiten des Gases und des Kalibriergases unter Normbedingungen abhangigen 50 
Korrekturfunktion bestimmt wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Korrek- 
turfunktion bestimmt wird: 

Pkorr = Pbda ■ [1 + K Pn • Ap n ], 55 
wobei Kp n eine Konstante oder eine Funktion von Pbda und/oder der Temperatur ist, und 

Apn = Pn,mess - Pnjcal- 



60 



65 
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4. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB auBerdem der Brennwert des Gases und des Kalibrier- 
gases unter Normbedingungen (Hn, me ss> H n ,kai) erfaBt wird und daB die Korrekturfunktion, die zur Bestimmung der 
unter Betriebsbedingungen gemessenen Dichte verwendet wird, vom Brennwert des Gases und des Kalibriergases 

5 unter Normbedingungen abhangig ist. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Korrek- 
turfunktion bestimmt wird: 

Pkorr = PBDA ' [1 + Kpn * A Pn + K H a " AH n ], 

10 

wobei Kp n und K Hn Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, und 

Ap n = Pn,mess ~ Pnjcal Und 
AH n ~ H njness - Hn ^al- 

15 

6. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Korrek- 
turfunktion bestimmt wird: 

Pkorr = Pbda • [1 + Kc n • Ac n + K Hn • AH n ], 

20 

wobei Kc n und K Hn Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, und 
Ac n = c n>racss - c fltkal und 

AH n = H n)rness - H n> kal- 

25 

7. Verfahren nach Anspruch 1 oder 4, dadurch gekennzeichnet, daB auBerdem die Temperatur des Gases unter Be- 
triebsbedingungen (T mess ) und die Temperatur des Kalibriergases beim Kalibrieren (T ka i) erfaBt werden und daB die 
Korrekturfunktion, die zur Bestimmung der unter Betriebsbedingungen gemessenen Dichte verwendet wird, von 
der Betriebstemperatur des Gases und der Kalibriertemperatur des Kalibriergases abhangig ist. 

30 8. Verfahren nach Anspruch 7, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Korrek- 

turfunktion bestimmt wird: 

Pkorr = Pbda • [ 1 + • Ap n + K Hn • AH n ] + K T • AT, 

35 wobei Kp n und K Hn Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, K T eine Konstante oder 

eine Funktion von Pbda ist, und 

Apn = pn,mess ~ Pnjcab 
AH n - Harness ~ H n> kal und 

40 AT = T raess -T kal . 

9. Verfahren nach Anspruch 7, dadurch gekennzeichnete daB die Dichte des Gases mit Hilfe der folgenden Korrek- 
turfunktion bestimmt wird: 

45 p korr = Pbda • [1 + • Ac n + K Hn • AH n ] + K T • AT, 

wobei K^ und Kh q Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, Kr eine Konstante oder 
eine Funkdon von Pbda ist, und 

50 Ac n = C nrness — C n Jkal> 

AH n = Harness - H„,kal und 

AT = T mess - Tkai. 

10. Verfahren nach einem der Anspriiche 1, 4 oder 7, dadurch gekennzeichnet, daB auBerdem der Kohlenstoffdi- 
55 oxidanteil des Gases und des Kalibriergases (xco2,mess> *co2,kai) erfaBt wird und daB die Korrekturfunktion, die zur 

Bestimmung der unter Betriebsbedingungen gemessenen Dichte verwendet wird, von dem Kohlenstoffdioxidanteil 
des Gases und des Kalibriergases abhangig ist. 

11. Verfahren nach Anspruch 10, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Kor- 
rekturfunktion bestimmt wird: 

60 

Pkorr = Pbda • [1 + Kp^ ♦ Ap n + K Hn • AH n + Kco2 ■ Ax C0 2] + K T • AT, 

wobei Kp n , Kh q und Kco2 Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, Kt eine Konstante 
oder eine Funkdon von Pbda ist, und 

65 

Ap n = Pn,raess ~ Pnjcab 
AH n = H njD css - Hn,kal> 

AT = T racss - T kal und 
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Ax C 02 = x C02.raess ~ *C02,kal- 

12. Verfahren nach Anspruch 10, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Kor- 
rekturfunktion bestimmt wird: 

Pkorr = Pbda ■ [1 + Kc n ■ Ac a + K Hn • AH n + Kc 02 • Axc 0 2] + K T • AT, 

wobei Kc n , Kh q und Kco2 Konstanten oder Funktionen von Pbda und/oder der Temperatur sind, Kt eine Konstante 
oder eine Funktion von Pbda ist, un d 

Ac n = C^mess — C^ai, 
AH n = Harness — H^fcai, 

AT = T mess - T kal und 

Ax c0 2 = X C 02,mess ~ *C02,kaL- 



Pkorr ~* Pbda 





+ K n 




X 








Pbda 




+ K , 




X 




1 + 1 ^ 

X 3 • c 

mesa 




Pkorr 



[( K io + P B0A • K T1 ) • (l^ ss • T^J] , 



BDA 



1 + 



■ C -l ^ C l*J " V PBOA I V llCal J. 



[(Kto + Pbda ' K n ) • (^ ess - T^J] 
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15 



13. Verfahren nach einem der Anspriiche 8 bis 12, dadurch gekennzeichnet, daB K T mit Hilfe der folgenden Funk- 
tion bestimmt wird: K T = Kt 0 + Pbda * K Tl , wobei K To und K Tl Konstanten sind. 

14. Verfahren zur Bestimmung der Dichte eines Gases, insbesondere von Erdgas, unter Betriebsbedingungen, wo- 
bei 

ein einen Schwingkorper enthaltender Betriebsdichteaufnehmer mit Hilfe eines Kalibriergases geeicht wird, 20 
und ein DichtemeBwert des Gases unter Betriebsbedingungen (Pbda) mi* Hilfe des Betriebsdichteaufhehmers ge- 
messen wird, dadurch gekennzeichnet, daB 

die Schallgeschwindigkeit und die Viskositat sowohl des Gases als auch des Kalibriergases jeweils unter Betriebs- 
bedingungen (c raess , c kal , Ti mess , Tl kal ) bestimmt wird und 

die Dichte des Gases unter Betriebsbedingungen aus dem DichtemeBwert mit Hilfe einer von der Schallgeschwin- 25 
digkeit und der Viskositat des Gases und des Kalibriergases abhangigen Korrekturfunktion bestimmt wird. 

15. Verfahren nach Anspruch 14, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Kor- 
rekturfunktion bestimmt wird: 



30 



35 



wobei T: Periodendauer des Schwingkorpers, 

T: Temperatur und 

Kc, K^, Kto, Kti Konstanten sind. 

16. Verfahren nach Anspruch 14, dadurch gekennzeichnet, daB die Dichte des Gases mit Hilfe der folgenden Kor- 40 
rekturfunktion bestimmt wird: 



45 



50 



wobei X: Periodendauer des Schwingkorpers, 

T: Temperatur und 

Kc, Kjj, Kto, Kti Konstanten sind. 

17. Verfahren nach einem der Anspriiche 1 bis 16, dadurch gekennzeichnet, daB Methan als Kalibriergas verwendet 
wird. 55 

18. Verfahren nach Anspruch 17, dadurch gekennzeichnet, daB Kc 2 und K n durch Kalibrierung mit einem zweiten 
Kalibriergas neben Methan, z. B. Argon, Ethen oder Stickstoff, bestimmt werden und als typenspezifische Konstan- 
ten behandelt werden. 

19. Verfahren nach einem der Anspriiche 14 bis 18, dadurch gekennzeichnet, daB die Schallgeschwindigkeit des 
Gases und/oder des Kalibriergases als Funktion der Dichte unter Betriebsbedingungen und der Dichte unter Norm- 60 
bedingungen und/oder der Temperatur und/oder des Brennwertes und/oder des Kohlenstoffdioxidanteils bestimmt 
wird. 

20. Verfahren nach Anspruch 19, dadurch gekennzeichnet, daB bei Bestimmung der Dichte von Erdgas als Kali- 
briergas Methan verwendet wird, 

daB die Schallgeschwindigkeit von Methan (CCH4) nach folgender Gleichung berechnet wird: 65 
CCH4 = A 0 + (A t + A 2 • T racss ) • T mcss + (A 3 + A4 • Pbda + A 5 • P 2 B da) ■ Pbda + + A 7 • T racss ) • T raess • Pbda 
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und die Schallgeschwindigkeit des Erdgases (c raess ) nach folgender Gleichung bestimmt wird: 

C meS s = Q:h4 + [(Bi + B 2 • Ap n + B 3 • T raess ) + (B 4 + B 5 - T mcss + B 6 • p BDA + B 7 • Pbda * Apn) ■ Pbda] ■ Ap n 
+ [(Ci + C 2 • T^) + (C 3 + C4 • Pbda + C 5 ■ T raess ) • p BDA ] • AH n + [D x • pbda] • X C o 2 , 

5 

wobei A0-A7, B L -B7, C1-C5 und Dt Konstanten sind, 
Xco 2 den Kohlenstoffdioxidanteil im Erdgas angibt, 

Ap n ~ p^nicss ~ Pn.CH* 

AH n = H n ,mess - H n> cH4- 

10 21. Verfahren nach einem der Anspruche 14 bis 20, dadurch gekennzeichnet, daB die Viskositat des Gases und/oder 

des Kalibriergases als Funktion der Dichte unter Betriebsbedingungen und der Dichte unter Normbedingungen und/ 
oder der Temperatur und/oder des Brennwertes bestimmt wird. 

22. Verfahren nach Anspruch 21, dadurch gekennzeichnet, daB bei Bestimmung der Dichte von Erdgas als Kali- 
briergas Methan verwendet wird, 
15 daB die Viskositat von Methan (1ICH4) nacn folgender Gleichung berechnet wird: 

TICH4 = Ao + A t • T + (A 2 + A 3 • Pbda) * Pbda 

und die Viskositat des Erdgases Cn me ss) nach folgender Gleichung bestimmt wird: 
Hmcss = ncm + Bi (1 - Pbda) • Ap n + B 2 • AH n 



20 



wobei A 0 -A 3 , Bi und B 2 Konstanten sind, 

Ap n = Pn,mess ~ Ph.CHa "nd 

25 AH n = H n mcss - H n>CH4 . 

23. Verfahren nach einem der Anspruche 1 bis 22, dadurch gekennzeichnet, daB zum Auffinden der Parameter und 
Konstanten geeigneter Funktionen, insbesondere der Korrekturf unktionen, mehrere MeBzyklen mit Referenzgasen 
durchgefuhrt werden, deren Zusammensetzungen bekannt sind. 



35 
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